Low-velocity anisotropic Dirac fermions on the side surface of topological insulators 



O 
(N 



Chang- Youn Moon, Jinhee Han, Hyungjun Lee, and Hyoung Joon ChoE] 
Department of Physics and IPAP, Yonsei University, Seoul 120-749, Korea 
(Dated: January 11, 2011) 

We report anisotropic Dirac-cone surface bands on a side-surface geometry of the topological 
insulator Bi2Se3 revealed by first-principles density-functional calculations. We find that the electron 
velocity in the side-surface Dirac cone is anisotropically reduced from that in the (lll)-surface Dirac 
cone, and the velocity is not in parallel with the wave vector k except for k in high-symmetry 
directions. The size of the electron spin depends on the direction of k due to anisotropic variation 
of the noncollinearity of the electron state. Low-energy effective Hamiltonian is proposed for side- 
surface Dirac fermions, and its implications are presented including refractive transport phenomena 
occurring at the edges of tological insulators where different surfaces meet. 
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Topological insulators (TIs) are characterized by odd 
number of Dirac-cone-like surface bands which guaran- 
tees robust metalhcity against arbitrary perturbations 
preserving time- reversal symmetry [l|, and by the 
helical spin texture of the surface states which leads 
to remarkable transport properties such as the absence 
of backscattering by non-magnetic defects as confirmed 
by experiments ^rm- The latter feature reminds us 
of graphene where charge carriers are also described as 
massless Dirac fermions, as demonstrated intriguingly by 
the Klein tunneling Q where normal incident wave ex- 
hibits perfect transmission through any finite potential 
barrier. A crucial difference is that TI involves real spins 
whereas graphene does pseudo-spins, so that, for exam- 
ple, one-dimensional wires of TI have maximum magne- 
toconductance at half integral flux quanta due to tt Berry 
phase acquired when circling around the wires 0, 0] ■ 

Promising TI materials [lo| include series of rhom- 
bohedral materials Bi2Se3, Bi2Te3, and Sb2Te3 proposed 
by a recent first-principles calculation with the ex- 
istence of single Dirac cone on their (111) surfaces in- 
deed confirmed by angle-resolved photoemission spectro- 
scopies ll|,ll2[. Among them, Bi2Se3 is considered most 
promising in the application point of view due to its rel- 
atively large bulk band gap of ^300 meV that can evade 
bulk conduction of thermally excited charge carriers at 
room temperature. This material has a layered struc- 
ture where the Bi2Se3 quintuple layers (QLs) are stacked 
along the rhombohedral (111) direction [13| combined by 
van der Waals type interaction. Thus, a natural cleav- 
ing face is the (111) surface between two QLs, and it has 
been the only surface considered theoretically or experi- 
mentally so far 14 1 . 

Surfaces other than the (111) surface, and interfaces 
between different surfaces, appear inevitably as step 
edges on cleaved (111) surfaces 0, lEI. Jisl a nd as side 
surfaces in nanowires or nanoribbons |16l4l9l |. Accord- 
ing to the topological band theory 0, [1] , a strong TI like 
Bi2Se3 should have odd number of Dirac cones on any 
arbitrary surface. However, existence of Dirac cones on 




FIG. 1: (Color online). Atomic structure of the Bi2Se3 (221) 
surface, (a) Side view with perspective, showing slanted QLs. 
(b) Top view of the topmost atomic layer, with surface unit 
vectors denoted by ai and a2. Using the experimental bulk 
structure the length of ai is 4.14 A, the j/-component of 
a2 is 11.26 A long, and the QLs are slanted by 58°. 



a surface other than the (111) surface and their physical 
properties have never been studied by realistic electronic- 
structure calculations or by direct measurements. 

In this Letter, we investigate the electronic structure 
of a side surface of Bi2Se3 using first-principles density- 
functional calculations, and find a Dirac-cone band struc- 
ture, for the first time on a surface other than the (111) 
surface, consistently with the Z2 topological band theory. 
We find that the side-surface Dirac cone has lower sym- 
metry, with elliptic equi-energy contour lines, and smaller 
velocity, introducing the refractive index in analogy with 
optics. We also present the low-energy Hamiltonian for 
the side-surface Dirac cone, and reflection and transmis- 
sion properties at interfaces between different surfaces. 

Our first-principles calculations are based on the 
density-functional theory with the generalized gradient 







ilBiE 







A FY 



FIG. 2: (Color online). Electronic structure of the Bi2Se3 
(221) surface, (a) Calculated side-surface electronic band 
structure from the BZ center to BZ boundaries along the two 
reciprocal lattice vectors. Filled regions denote bulk contin- 
uum states projected on the surface BZ. (b) Side view and 
(c) top view of a wavefunction at the Dirac point, which show 
that the state is localized near the surface. 

approximation [20| and ab-initio norm-conservinff pseu- 
dopotentials as implemented in SIESTA code [21[. We 
employ partial-core correction for Bi and Se pseudopo- 
tentials, pseudo-atomic orbitals (double-^-polarization 
basis set) for electronic wavefunctions, and a cutoff en- 
ergy of 500 Ry for real-space mesh. We use a supercell 
containing a slab of Bi2Se3 and 9 A thick vacuum to 
study surface states while preventing interaction between 
periodic images. Brillouin-zone (BZ) integration is per- 
formed with a two-dimensional (2D) 12 x 6 A:-point grid. 
Spin-orbit interaction is treated by relativistic pseudopo- 
tentials in the fully nonlocal form 22]. Atomic positions 
are taken from the experimental bulk structure [ISj . 

Figure 1 shows the atomic structure of the side sur- 
face considered in our present work. This side surface is 
stoichiometric in the sense that the topmost atomic layer 
of the surface and other atomic layers in parallel with 
it have stoichiometric Bi and Se compositions for Bi2Se3 
within each atomic layer. The normal direction of the 
side surface is (221) in Miller indices of the rhombohedral 
bulk structure, making an angle of 58° to the (111) plane. 
In Fig. 1, no bond is drawn between Se atoms belonging 
to different QLs to represent the (111) cleaving planes 
with real-space gaps which are along the x direction and 
repeated in the y direction. The (221) surface has lower 
symmetry than the (111) surface, exhibiting anisotropy 
between x and y directions clearly. We construct a slab 
of eleven atomic layers of Bi2Se3 with the slab thickness 
of ~^30 A to study surface electronic structures. 

Figure 2(a) shows calculated electronic band structure 
along two lines in 2D BZ: from F to A and from F to Y. 
Here, A and Y arc BZ boundaries along the reciprocal 
lattice vectors bi and b2, respectively, associated with 
the real-space unit vectors ai and a2 shown in Fig. 1(b). 
In Fig. 2(a), a crossing of linear bands is clearly visible at 
the r point with the Dirac-point energy coinciding with 
the Fermi energy (Ef). The bulk band continuum is pro- 
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FIG. 3: (Color online). The Dirac cone on the Bi2Se3 (221) 
surface, (a) Calculated 2D electronic band structure near the 
Dirac point at k = 0. Equi-energy contours for -150, -100, 
-50, 50, 100, and 150 meV from the Dirac point are drawn on 
the energy surface and then projected onto the bottom plane, 
(b) Equi-energy contour at 20 meV below the Dirac point. 
Arrows denote calculated spins of the states at each k. 

jected onto the 2D surface BZ as filled areas in Fig. 2(a), 
where the bulk band gap is about 300 meV, consistently 
with previous calculation Q. To check if the crossing 
bands are indeed surface states, we plot the squared am- 
plitude of the wavefunction, ['(/'(r)^, for a state at the F 
point in Figs. 2(b) and (c). We find that it is mainly 
located on the surface atoms and decays rapidly into 
the bulk region, showing the surface-state nature of the 
crossing bands. Around the F point in Fig. 2(a), we no- 
tice that the slope of the linear band is steeper, i.e. the 
Fermi velocity is larger, along the T-A direction than F- 
Y, indicating anisotropy of the Dirac cone. Particle-hole 
asymmetry is also visible, where band linearity persists 
down to ~200 meV below Ef while the two bands bend 
significantly only ^30 meV above Ef. 

The shape of the side-surface Dirac cone is more obvi- 
ous in the 2D plot of the bands near the F point, as shown 
in Fig. 3(a). The Dirac cone is anisotropic, with the slope 
of the band and hence the Fermi velocity being larger in 
the X direction than in the y direction. As mentioned 
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above, the linear band feature holds only within a very 
small energy window near Ef, so the energy contours of 
50 meV above Ef (in black ) and below Ef (in orange) 
look very different from each other [Fig. 3(a)]. At ener- 
gies near Ef close enough to ensure the band linearity, 
the energy contours are basically very close to ellipses. 

We obtain the spin of a Dirac-cone state on the (221) 
surface by evaluating expectation values of the three 
spin-operator components, and the results are shown in 
Fig. 3(b). As expected for the helical spin fermion, the 
spin is almost perpendicular to the wave vector k. Mean- 
while, since each electronic state in this material is not a 
spin eigenstate but a noncollinear spin state due to the 
spin-orbit interaction, the spin size is not necessarily 1 
in general [i^. Unlike the (111) surface where the 
spin size is independent of k, the spin size in the (221)- 
surface Dirac cone is found to depend on k, having its 
maximal value of 0.83 /is at = and gradually de- 
creasing toward 0.55/isatfcj^=0 around the Dirac cone 
[Fig. 3(b)]. 

Now we consider low-energy effective Hamiltonian sim- 
ilar to the previously suggested one ^ which basically 
corresponds to a Rashba-type surface spin-orbit Hamil- 
tonian. Considering the reduced symmetry, the (221)- 
surface Hamiltonian in the first order of k would be 



(1) 



where ax and ay are the Pauli matrices. The Fermi veloc- 
ities Vx and Vy in x and y directions on the (221) surface 
are 3.1 x 10^ m/s and 1.4 x 10^ m/s, respectively, from 
our first-principles calculations. Energy eigenvalues are 



E^±hjvlkl 



■v^k^, 



and the corresponding eigenstates are 

„ik r / 1 



^k(r) 



^ y±{-iVxkx + Vyky)/^Jvlkf+vp{] 



(2) 



(3) 



In Eqs. (2) and (3), the plus sign is for the conduc- 
tion band while the minus sign is for the valence band. 
The equi-energy contour from Eq. (2) is an ellipse, which 
agrees with Fig. 3(b) a.t E ~ —20 meV. 

As mentioned earlier, interfaces of different surfaces 
are very often involved in many experimental situations 
where two different surfaces meet as a form of step 
edges on a cleaved TI surface or geometric edges in TI 
nanowires and nanoribbons. As our result shows that 
different TI surfaces have Dirac cones with different ve- 
locities, we can expect reflection and refraction behaviors 
in electronic transports at the interfaces of TI surfaces in 
close analogy with optics with distinct refractive-index 
materials. The concept of the refractive index can be in- 
troduced to a TI surface, and it provides a new insight 
for electronic transport in various TI geometries. We dis- 
cuss two basic simple interface of TI surfaces and 
a Fabry-Perot interferometer. 
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FIG. 4: (Color online). Ray optics of Dirac fermions in TI. 
(a) Schematic diagram of the edge of the (111) and (221) 
surfaces in Bi2Se3. (b) Transmission probability T and re- 
fraction angle d as functions of incident angle </!> in (b). (c) 
Schematic diagram of the Fabry-Perot interferometer consist- 
ing of (111)/(221)/(111) surfaces, (d) T as a function of </> 
with the (221)-surface widths of 10, 50, and 100 nm in (c). 



First we consider a simple interface between (111) 
and (221) surfaces which is parallel to the x direction 
[Fig. 4(a)]. When a Dirac-cone state of the wave vector 
k is incident from the (111) surface to the interface at 
?/ = 0, with the incident angle 0, we can obtain the trans- 
mission probability T and the refraction angle d from 
the equality of x-components of incident, reflected, and 
transmitted wave vectors, the energy conservation, and 
the continuity of the wave functions at the interface, sim- 
ilarly to the graphene case previously studied [1, [2J, [2^ : 



9') + cos(0 -f 0') 



1 + cos{(t) + 9') 



(4) 



where 9' = tan~^ (vxkx/vyqy) appears from the rela- 
tive phase of the two spinor components in Eq. (3). 
The refraction angle 9 is given hy 9 = tAn~^ (kx/qy), 
where k^ = k simp and qy = v'^k'^ — v'^k'^ / Vy is the 
y-component of the wave vector in the (221) surface re- 
gion. The Fermi velocity vq in the (111) surface is about 
5 X 10^ m/s in our calculation, consistently with previous 
work [§]. Here we note that the incident and refraction 
angles satisfy the Snell's law with the refractive index n 
of the TI surfaces defined as the inverse of the phase ve- 
locity {vph = E/hk) of Dirac fermions, which, in general, 
depends on the k direction but not on the size of k. 

Calculated transmission probability T and refraction 
angle 9 are plotted in Fig. 4(b) as functions of incident 
angle (p ai E — 2Q meV. At normal incidence (0 = 0), T 
is 1, which means a perfect transmission, demonstrating 
the absence of back scattering due to the spin texture. 
As (f) increases to 90°, T slowly decreases to 0, being 
greater than 0.5 for (f) less than 80°. At the same time, the 
refraction angle 9 increases slowly from 0° and reaches 
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20° at (j) = 90°. Thus, in the opposite situation in which 
the Dirac fermion is incident from the side surface to the 
(111) surface, the refraction angle reaches 90° already at 
the incident angle of 20° and at larger incident angles 
the Dirac fermion cannot transmit through the interface. 
This implies that the total internal reflection can occur 
when the electron is incident from the side surface to the 
(111) surface, similarly to the optical system where the 
incident wave resides in a dielectric with higher refractive 
index than the other side of the interface. This opens a 
chance to produce evanescent waves of Dirac fermions on 
the (111) surface and their related phenomena. 

Another relevant geometry is a Fabry-Perot interfer- 
ometer where the (221) surface is made between two 
(111) surfaces [Fig. 4(c)], which might be a possible can- 
didate to describe the step edges on the (111) surface in 
many previous experiments [j, Q . With boundary con- 
ditions similar to the above interface problem, we obtain 
reflectivity r and hence the reflectance R — |rp in anal- 
ogy with the graphene case Q 

^ _ sm^{qyD){sm(p - s'md')'^ 

cos^i/) cos^O' + sm^{qyD){sm<l) — smO')"^ ' 

where D is the width of the (221) surface region 
[Fig. 4(c)], and the other parameters are defined simi- 
larly to the simple-interface case above. 

The transmission probability through the interferome- 
ter, T = 1 — i?, is plotted as a function of the incident an- 
gle (j) in Fig. 4(d), considering three different side-surface 
widths D = 10, 50, and 100 nm, with £: = 20 meV. For 
all cases, we have the perfect transmission for the nor- 
mal incidence ((/> = 0) as in the simple-interface case, 
because of the forbidden back scattering. We also expect 
a perfect transmission to occur resonantly when the de- 
structive interference happens between the incident and 
reflected waves, i.e., QyD — rnr for positive integer n from 
Eq. (5). Because of relatively low Fermi velocity in the 
side surface region, _D = 10 nm is not large enough for 
the resonant transmission to occur at i? = 20 meV, so 
T monotonically decreases from 1 to with increasing 0. 
For larger widths, 50 and 100 nm, T reaches 1 at = 53°, 
where the resonance condition is fulfilled. In the case of 
D = 100 nm, the resonance at 4> = 53° is narrower, with 
less transmission at (f> between 0° and 53°. 

Our Fabry-Perot structure can also be considered as an 
electronic analogy of the 'optical fiber' where propagating 
waves are confined within the high refractive index (low 
Fermi velocity) region by the total internal reflection. If 
the side surface region is clean enough to prevent the 
waves from escaping into the (111) surface region with 
high incident angles by defect scattering, this structure 
might be used as a one-dimensional conducting wire. 

In conclusion, we performed first-principles calcula- 
tions to investigate the electronic structure of a side sur- 
face of topological insulator Bi2Se3, finding a Dirac cone 



on the side surface, consistently with the topological- 
band-theory argument. The side-surface Dirac cone is 
found to be slow and anisotropic in an elliptic shape, 
with the spin size dependent on k. The low-energy effec- 
tive Hamiltonian is proposed for the side-surface Dirac 
fermions, and the concept of the refractive index is intro- 
duced to the TI surface. Refraction, reflection, and trans- 
mission behaviors of Dirac fermions are demonstrated in 
the simple TI interface and the Fabry-Perot TI interfer- 
ometer. Our idea of making TI junctions using different 
surfaces can be applied to form more complex structures 
such as superlattices, opening a chance for novel physical 
properties of Dirac fermions. 
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